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1 How it all began 

I started collecting integer sequences in December 1963, when I was a graduate 
student at Cornell University, working on perceptrons (or what are now called 
neural networks). Many graph-theoretic questions had arisen, one of the simplest 
of which was the following. 

Choose one of the n n_1 rooted labeled trees with n nodes at random, and 
pick a random node: what is its expected height above the root? To get an 
integer sequence, let a n be the sum of the heights of all nodes in all trees, and 
let W n — a n /n. The first few values W\, W2, ■ ■ ■ are 

0, 1, 8, 78, 944, 13800, 237432, 



a sequence engraved on my memory. I was able to calculate about ten terms, but 
I needed to know how W n grew in comparison with n™, and it was impossible 
to guess this from so few terms. So instead I tried to guess a formula for the 
nth term. Again I was unsuccessful. Nor could I find this sequence in Riordan's 
book Q, although there were many sequences that somewhat resembled it. 

So I started collecting all the sequences I could find, entering them on 
punched cards, thinking that if any of these sequences came up in another prob- 
lem, at least I would know what they were. 

I never did find that sequence in the literature, but I learned Polya's theory 
of counting and (with John Riordan's help) obtained the answer, which appears 



in |87[ ]. There is a simple formula, although maybe not simple enough to be 
guessed: 

n-2 k 

^n = (»-l)!EF- 

fc=0 

An old formula of Ramanujan | p9| then implies that 



as n 



which was what I needed. That sequence became number A435 in the collection. 

The idea of a "dictionary" of integer sequences was received with enthusiasm 
by many people, and in 1973 I published containing about 2400 sequences, 
arranged lexicographically. One correspondent commented on the book by saying 
"There's the Old Testament, the New Testament and the Handbook of Integer 
Sequences" . 

Over the next twenty years an enormous amount of new material arrived, 
preprints, reprints, postcards, typewritten letters, handwritten letters, etc., and 
it was not until 1995 that — with Simon Plouffe's help — a sequel |)3| appeared. 
This contained 5500 sequences. 

Around the same time I set up two services that can be used to consult the 



database via electronic mail. The first has the address jequences@research.att.com , 



and simply looks up a sequence in the table. The second email address, which 



is 3uperseeker@research.att.com, tries hard to find an explanation even for a 
sequence not in the table. 

A large number of new sequences started arriving as soon as Q appeared, 
and when in 1996 the total number reached 16,000, three times the number in 
the book, I decided to set up a web site for the database Q. New sequences 
still continue to pour in, at about 10000 per year. At present, in April 1999, the 
database contains about 50000 sequences, and the web site receives over 2500 
hits per day. 

The main reason for this rapid expansion is that in the two books I only 
included sequences that had been or were about to be published. For the on- 
line version, where storage space is no longer a limitation, any well-defined and 
sufficiently interesting sequence is eligible for inclusion. 

There is now also an electronic Journal of Integer Sequences flc| and a mailing 
list for sequence fans |89) . 

The following sections describe how the database is used (Section 2) and 
the kinds of sequences it contains (Section 3). Section 4 discusses a few "hard" 
sequences and Section 5 some recursive examples. Then Sections 6-8 describe 
sequences associated with meandering rivers and stamp-folding, extremal codes 
and lattices, and Levine's sequence. 

Besides integer sequences, the table also contains many examples of arrays 
of numbers, Pascal's triangle, for example — see Section 9. The final three sec- 
tions discuss sequences associated with the Wythoff array, the boustrophedon 
transformation of sequences, and Tchoukaillon solitaire. 
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2 How the database is used 



The main applications of the database are in identifying sequences or in finding 
out the current status of a known sequence. 

The database has been called a mathematical analogue of a "fingerprint 
file" JnJ. You encounter a number sequence, and wish to know if anyone has 
ever come across it before. If your sequence is in the database, the reply will 
provide a description, the first 50 or so terms (usually enough to fill three lines 
on the screen), and, when available, formulae, recurrences, generating functions, 
references, computer code for producing the sequence, links to relevant web sites, 
etc. 

Let me illustrate how the database is used with a typical story. Last summer 
the following question arose at AT&T Labs in connection with a quantization 
problem ^Kj. Given an n-dimensional lattice A, for which integers TV does A 
have a sublattice of index TV that is geometrically similar to TV? 

For the two-dimensional root lattice A 2 , for example, it is easy to see that a 
necessary and sufficient condition is that TV be of the form a 2 + ab + 6 2 . 

However, for the four-dimensional lattice A4 the situation is more compli- 
cated. The first thing we did was to run a computer search, which showed that 
A4 has a similar sublattice of index TV = c 2 if and only if c is one of the numbers 



1, 4, 5, 9, 11, 16, 19, 20, 25, 29, 31, 36, . 



This turned out to be sequence A31363 in the database, with a reference to 
Baakc (3), where it appears in an apparently different context as the indices 
of coincidence site sublatticcs in a certain three-dimensional quasicrystal. 
identifies these numbers as those positive integers in which all primes congruent 
to 2 or 3 (modulo 5) appear to an even power. This was a very useful hint in 
getting started on our problem. We were able to show that this is the correct 
condition for the A4 lattice, and to find analogous results for a number of other 
lattices |l5). (However, we have not yet found a direct connection between A4 
and the quasicrystal problem. Nevertheless, the occurrence of the same numbers 
in the two problems cannot be entirely coincidental.) 

My files contain many similar stories. Some other examples can be found in 
Chapter 3 of @. 

Most of the applications however are less dramatic. One encounters a se- 



quence in the middle of a calculation, perhaps 



1,1,2,3,6,10,20,35,70,126, 252, . . ., 

and one wants to know quickly what it is — preferably a formula (in this case 
it is ([„™ 2 ])) or g enera ting function. 

It is worth emphasizing a special case of this: the simplification of binomial 
coefficient sums. Powerful methods are available for simplifying such sums by 
computer J7§, @- But if one is in a hurry, one can first try evaluating the 
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initial terms and looking up the sequence in the table. E.g. the sum 

, . (In - 2k\ 2 /2k 

produces 

1,8,88,1088,14296,195008,... 
( A3G917] ), and the entry in the table supplies a recurrence 

n 3 a(n ) = 16 (n - (2n 2 - In + l)a(n - 1) - 256(?i. - l) 3 a(n - 2) 

and a reference (fn\, p. ix, although the sum is misstated in the first printing). 
I have begun entering into the database all the sequences corresponding to the 
singly-indexed identities in Gould's table |32| of binomial coefficient identities. 

A related application is in identifying arithmetic inequalities. Suppose you 
suspect that 

a(n) > d(n) + 0(n), for n > 2 , (1) 

where a, d and 4» are respectively the sum of divisors, number of divisors, and 
Euler totient functions. You could evaluate the sequence formed by the difference 
of the two sides, which for n > 1 is 

-1, 0, 0, 2, 0, 6, 0, 7, 4, 10, 0, 18, 0, 14, 12, . . . 

( |A46520D . The table then points you to a reference ( ]73|| , §1.3.1) where this is 
stated as a theorem. (Again I would like to get more examples of such sequences.) 

Another important use for the database is in finding out the present state 
of knowledge about some problem. A second story will illustr ate th is. The num- 
ber of Latin squares of order n (Q n ) is given by sequence A31£ (see Fig. |l|). 
Computing Q n is one of the famous hard problems in combinatorics. 



1 1 

2 1 

3 1 

4 4 

5 56 

6 9408 

7 16942080 

8 535281401856 

9 377597570964258816 
10 7580721483160132811489280 



Fig. 1. Sequence A315 , the number of Latin squares of order n (McKay and 
Rogoyski jn|). 
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In 1991 Brendan McKay, at the Australian National University in Camberra, 
computed Qio (which as can be seen is a rather large number). When he checked 
the database he found that the same value had recently been obtained by Eric 
Rogoyski of Cadence Design Systems in San Jose, California. As it turned out 
the two methods were similar but not identical, and he and Rogoyski ended 
up writing a joint paper JnJ, acknowledging the database for bringing them 
together. 

3 Types of Sequences 

The database contains sequences from all branches of science, including 

— enumeration problems (combinatorics, graph theory, lattices, etc.) 

— number theory (number of solutions to x 2 + y 2 + z 2 = n, etc.) 

— game theory (winning positions, etc.) 

— physics (paths on lattices, etc.) 

— chemistry (sizes of clusters of atoms, etc.) 

— computer science (number of steps to sort n things, etc.) 

— communications (m-sequences, weight distributions of codes, etc.) 

— puzzles 

— etc. 

To be accepted into the database, a sequence must be well-defined and inter- 
esting. However, when in doubt, my tendency is to accept rather than to reject. 
The amazing coincidences of the Monstrous Moonshine investigations make 
it difficult to say that a particular sequence, no matter how obscure, will never 
be of interest. 

Readers are urged to send me any sequences they come across that are not 
at present in the database. There is a convenient electronic form for this purpose 
in f9^| . Some of the reasons for sending in your sequence are as follows: 

— this stakes your claim to it 

— your name is immortalized 

— the next person who comes across it will be grateful and, not least, 

— you may benefit from this yourself, when you come across the same sequence 
some weeks from now. 

Often one finds that a particular project may involve dozens of sequences, all 
variants of a few basic ones. Ideally you should send them all to the database! 

4 Hard sequences 

One of the keywords used in the database is "hard", which indicates that the 
term following those given is not known. Besides the Latin squares problem men- 
tioned above, some other classic hard sequences are the following: 
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Projective planes. The number of projective planes of orders n = 2, 3, . . . , 10 
( |A1231| ): 

1,1,1,1,0,1,1,4,0, 

where the last term refers to the result of Clement Lam et al. | | (completing 
work begun in p^] ) that there is no projective plane of order 10. 

The Poincare conjecture. The number of differential structures on the n- 
sphere, for n = 1, 2, . . . , 16 ( |A1676| ): 

1, 1, 1?, 1, 1, 1, 28, 2, 8, 6, 992, 1, 3, 2, 16256, 2, 

as given by Kervaire and Milnor The Poincare conjecture is that the third 
term is 1. 



Dedekind's problem. The number of monotone Boolean functions of n vari- 
ables, for n = 1, 2, 3, . . . , 8 ( |A7153| ): 

1,4, 18, 166, 7579, 7828352, 2414682040996, 
56130437228687557907786 



where the last term was computed by Wiedemann JlOOj] . (As is the case for 
many of these examples, there are several other versions of this sequence in the 
database.) 



The Hadamard maximal determinant problem. What is the max imal dc - 
terminant of an n x n {0, l}-matrix? The values for n = 1, 2, . . . , 13 are ( |A3432| ): 

1, 1, 2, 3, 5, 9, 32, 56, 144, 320, 1458, 3645, 9477, 

where the last two terms are due to Ehlich, and Ehlich and Zeller [^5| , |2^| . For 
n = — 1 (mod 4), Hadamard showed that the nth term is equal to 

(n + l) ( " +1) / 2 /2" , 

provided that what is now called a Hadamard matrix of order n exists. In some 
cases conference matrices give the answer when n = 1 ( mod 4), but the problem 
of finding the other terms in the sequence seems to have been untouched for over 
35 years. It would be nice to have confirmation of the above values as well as 
some more terms! 



Enumerating Hadamard matrices. The number of Hadamard matrices of 
orders n = 4, 8, 12, . . . , 28 ( |A7299[ ) is 

1,1,1,5,3,60,487, 



where the last entry is the work of Kimura |5(J , |5lJ] , |^2j] , |53| . The Hadamard 
conjecture is that such a matrix always exists if n is a multiple of 4. Judging by 



6 



how rapidly these numbers are growing, this should not be hard to prove, yet 
it has remained an open question for a century. Of course, as the example in 
Section 7 shows, such numerical evidence can be misleading. 



The kissing number problem. How many spheres can touch another sphere 
of the same size? For arrangements that occur as part of a lattice packing, the 



answers are known for n = 1, 2, . . . , 9 ( A1116 ) 



2,6,12,24,40,72,126,240,272, 

the last term being due to Watson (see |l6|). For nonlattice packings, all we 
know is 

2, 6, 12, ?,?,?,?, 240, > 306 . 
The best bounds known in dimensions 4, 5, 6 and 7 are respectively 

24-25, 40-46, 72-82 and 126-140 

— see [flql for further information. 



Sorting by prefix reversal. If you can only reverse segments that include 
the initial term of the current permutation, how many reversals are needed to 
transform an arbitrary permutation of n letters to the identity permutation? To 
state this another way |24j : 

The chef in our place is sloppy, and when he prepares a stack of pancakes 
they come out all different sizes. Therefore, when I deliver them to a 
customer, on the way to the table I rearrange them (so that the smallest 
winds up on top, and so on, down to the largest at the bottom) by 
grabbing several from the top and flipping them over, repeating this 
(varying the number I flip) as many times as necessary. If there are n 
pancakes, what is the maximum number of flips (as a function f(n) of 
n) that I will ever have to use to rearrange them? 

The only exact values known are /(!),...,/ (9): 



0,1,3,4,5,7,8,9,10, 

due to Garey, Johnson and Lin, and Robbins (see Q, p9|] ). It is also known that 
f(n) > n+l forn > 6, f(n) > 17n/16 if n is a multiple of 16 (so /(32) > 34), and 
f(n) < (5n + 5)/3, the last two bounds^ being due to Gates and Papadimitriou 
p9| . Again it would be nice to have more terms. 

1 (Added May 5, 2000.) The bound f(n) < (5n + 5) /3 was independently obtained by 
E. Gyori and G. Turan, Stack of pancakes, Studia Sci. Math. Hungar., 13 (1978), 
133-137. I thank Laszlo Lovasz for pointing out this reference. 
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Note added Jan. 17, 2001. John J. Chew III (Department of Mathematics, 
University of Toronto) has found that f(!0) through /(13) are 11,13,14,15, 
respectively, so the sequence begins: 

0,1,3,4,5,7,8,9,10,11,13,14,15 

This has now been added to the database as sequence |A058986|. 



5 Recursive sequences 

Whereas the sequences in the previous section enumerated some class of objects, 
the following are self-generated. 



Differences = complement ( A5228| ): 

1, 3, 7, 12, 18, 26, 35, 45, 56, 69, 83, 98, 114, 



The differences ( |A30124| ) 

2, 4, 5, 6, 8, 9, 10, 11, 13, 14, 15, 16, 17, 19, . . . 

are the terms not in the sequence! This is one of many fine self-generating se- 
quences from Hofstadter |3q|. 



Golomb's sequence. The nth term is the number of times n appears (|Al462p : 

1,2,2,3,3,4,4,4,5,5,5,6,6,6,6,7, 7, 7,7,8,... 
The nth term is the nearest integer to (and converges to) 



where r = (1 + VE)/2 M, @ Section E25]. 



Wilson's primeth recurrence. a n+ i is the q n -t h pr ime ( A7097 ), shown in 
Fig. H The sequence was sent in by R. G . Wilson V [ 102 1 , and the last few terms 
were computed by P. Zimmermann [102] and M. Deleglise |Q. Their algorithm 
is a slightly speeded up version of an algorithm for computing ir(x), the number 
of primes not exceeding x, due to J. C. Lagarias, V. S. Miller and A. M. Odlyzko 
( see Q] ) . It is quite remarkable that it is possible to compute so many terms of 
this sequence. 

Recaman's sequences, (i) a n = a„_i — n if a„_i — n > and a n _i — n has 
not already occurred in the sequence, otherwise a n = a n _i + n ( A5132| ): 



1, 3, 6, 2, 7, 13, 20, 12, 21, 11, 22, 10, 23, 9, 



(ii) a„+i = a„/n if n divides a„, otherwise a„+i = na„ ( A8336): 
1, 1, 2, 6, 24, 120, 20, 140, 1120, 10080, . . . 
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1 

2 
3 
5 
11 
31 
127 
709 
5381 
52711 
648391 
9737333 
174440041 
3657500101 
88362852307 
2428095424619 
75063692618249 
2586559730396077 

Fig. 2. a rl +i is the a„-th prime. 



These were sent in by B. Recaman |84|| . How fast do they grow? 

The $10,000 sequence. In a colloqui um talk at AT&T Bell Labs Jlf], John 



Conway discussed the sequence (A4001) 

1,1,2,2,3,4,4,4,5,6,7,7,8,8,8,8,9,... 
defined by (for n > 3) 

a(n + 1) = a(a(n)) + a(n + 1 — a(n)) . 

(In words, a(n + 1) is the a(n)th term in from the left plus the a(n)th term in 
from the right.) This sequence seems to have been introduced by either David 
Newman or Douglas Hofstadter around 1986. In his talk, Conway said that he 
could prove that — ► A, and offered $10,000 for finding the exact n at which 
— \ last exceeds My colleague Colin Mallows did not take long to 

analyze the sequence, and came up with an answer of 6083008742 §7$. 

Colin tells me that in fact the problem is actually much easier than either he 
or John Conway had believed, and the true answer is 1489. A recent paper by 
Kubo and Vakil Q also studies this sequence and its generalizations. See also 
[H Section E31]. 

Many variants of this sequence have not yet been analyzed. Even the rate of 
growth of this one is not known: a(l) = a(2) = 1, a(n) = a(a(n — 2)) + a(n — 



a(n - 2)), A5229 M, which begins 



1, 1, 2, 3, 3, 4, 5, 6, 6, 7, 7, 8, 9, 10, 10, 10, 11, 12, 12, . 
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The Prague clock sequence. This is not really recursive in the same sense as 
the preceding sequences, but it seems to fit in here. It is, in any case, delightful. 



This is |A28354| : 1, 2, 3, 4, 32, 123, 43, 2123, 432, 1234, 32123, 43212, 34321, 23432, 
123432, 1234321, 2343212, 3432123, 4321234, 32123432, 123432123, 43212343, 
2123432123, 432123432, 1, 2, 3, 4, 32, . . . According to Q, the sequence indicates 
how the astronomical clock in Prague strikes the hours. There is a single bell 
which (i) makes from 1 to 4 strokes at a time, (ii) the number of strokes follows 
the sequence 

...3212343212343... , 

(hi) at the nth hour, for n = 1,2, ... , 24, the strokes add to n, and (iv) at the 
25th hour there is a single stroke (so the sequence has period 24). As the reader 
will see by studying the sequence, its existence depends on two coincidences! 

Cayley's mistake. Since the sequences in the database arc numbered Al, 
A2, A3, . . ., several people humorously proposed the "diagonal" sequence (now 
A31135| ) m which the nth term is equal to the nth term of Ar 



in: 

1, 2, 1, 0, 2, 3, 0, 6, 8, 4, 63, 1, 316, 42, 16, . 



and the even less well-defined sequence (now A37181) with nth term equal to 



A31135 (n) + 1. I resisted adding these sequences for a long time, partly out of a 
desire to maintain the dignity of the database, and partly because ^ was only 
known to 11 terms! 



Sequence A 2 gives the number of "centered hydrocarbons with n atoms", 
and is based on an 1875 paper of Cayley g|. The paper is extremely hard to 
follow, and gives incorrect values for n = 12 and 13. The errors it contains were 
reproduced in §, |»J and even though Herrmann |37| had pointed out 



these errors in 1880. As far as we can tell, a correct verion of this sequence was 
never published until Eric Rains and I did so in 1999 |83|]. Although Henze and 
Blair, Polya and many others have written articles enumerating related families 
of chemical compounds (see p3j for a brief survey), this sequence seems to have 
been forgotten for over 100 years. 

Once we determined what it was that Cayley was trying to enumerate, Polya's 
counting theory quickly gave the answer, and the correct version of 



A22 



is now 



in the database (as are the two diagonal sequences mentioned above - mostly 
because users of the database kept proposing them). 



6 Meanders and stamp-folding 

The meandric and stamp-folding numbers are similar to the better-known Cata- 
lan numbers (sequence A108) in that they are fundamental, easily described 
combinatorial quantities that arise in many different parts of mathematics, but 
differ from them in that there is no known formula, and in fact seem to be quite 
hard to compute. 
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The stamp-folding problem has a history going back to at least Lucas |i3| 
in the nineteenth century, while the meandric problem seems to have been first 
mentioned by Poincare ]80| ]. The meandric problem asks: in how many different 
ways can a river (starting in the South- West and flowing East) cross a road n 
times? For n = 5 crossings there are M5 = 8 possibilities, shown in Fig. |^. The 



sequence Mi, M 2 , M 3 , . . . ( A5316) begins 



1, 1, 2, 3, 8, 14, 42, 81, 262, 538, 1828, . . . 

These are called meandric numbers, since the river meanders across the road. 
The even- numbered terms M2, M±, Mq, . . . give the number of different ways an 





7^/ 















Fig. 3. The eight ways a river (going from South- West to North-East) can cross 
a road five times. 



oriented line can cross a Jordan curve ( A531£ ) . There are several other interpre- 
tations, one of which is the number of "simple alternating transit mazes" |7^| , 

s 

The stamp-folding problem asks the same question, but now the line is only 
semi-infinite. Equivalently, how many ways are there to fold a strip of n stamps? 
Pictures illustrating the first few terms of the stamp-folding sequence (A1011) 
can be found on the front cover of |nj and in Figure M4587 of p3l. The sequence 
begins 

1, 1, 2, 5, 14, 38, 120, 353, 1148, 3527, . . . 
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No polynomial-time algorithm is known for computing either sequence (on 
the other hand it is not known that such algorithms do not exist). The best 
algorithms known require on the order of nC n steps, where C„ = -^fj( 2 ") is 




the nth Catalan number (A108): these algorithms are due to Koehler 
the stamp-folding problem and to Knuth and Pratt |54j and Reeds for the 
meandric problem. They are too complicated to describe here. 

Using these algorithms, Stephane Legendre |6l) has extended the stamp- 
folding sequence to 26 terms and the meandric sequence to 25 terms. Knuth and 
Pratt j54j] have computed 17 terms of the M 2 , M 4 , M 6 , . . . subsequence. Lando 
and Zvonkin @ and Di Francesco, Golinclli and Guitter @, @, 
have also studied these sequences. 

The exact rate of growth of these sequences is not known. The best bounds 
for Mm presently known appear to be due to Reeds, Shepp and Mcllroy It 
is easy to see that M% n is submultiplicative, and that C n < Mm < G^, which 
implies that 

H = lim M^ n 

n— *oo 

exists and satisfies 4 < « < 16. In |8j| it is shown that 

8.8 < n < 13.01 . 

Besides the papers already mentioned, the meandric and stamp-folding prob- 
lems have recently been discussed by Arnold jl[], Di Francesc o []l9| , Harris fl(| , 
Lando and Zvonkin |59), ||(J, Lunnon Q, Sade [p8| , Smith |95| and Touchard 



7 Extremal codes and lattices 

Let C be a binary linear self-dual code of length n in which the weight of every 
codeword is a multiple of 4, and let 



Wc(x,y) = y"V 



n— wt(c) wt(c) 

cEC 



be its weight enumerator. Examples are the Hamming code of length 8, with 
weight enumerator 

/ = x 8 + Ux 4 y 4 + y 8 , 
and the Golay code of length 24, with weight enumerator 

g = x 24 + 759x 1 V + 2576a; 1 V 2 + 759x 8 y 16 + y 24 . 

A remarkable theorem of Gleason says that the weight enumerator of any 
such code C is a polynomial in / and g. (References for this section arc 
[ p5| , E.g. if C has length n = 72, its weight enumerator can be written as 

W C - a a f + ai fg + a 2 fg 2 + a 3 g 3 , 
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for rational numbers ao, . . . , 03. If we choose these numbers so that the minimal 
distance of this (hypothetical) code C is as large as possible, we find that 



W c = 1 + Cte 4 + Ox 8 + Ox 12 + 249849a; 16 + 18106704a; 20 + 462962955a; 24 



so that C would have minimal distance 16 (the coefficients form sequence |A18236| ) 



The coefficients in this "extremal" weight enumerator are all nonnegative inte- 
gers, but whether a code exists with this weight enumerator is an important 
unsolved question. 

One can perform this calculation for any length that is a multiple of 8, and in 
[|70l it was shown that when as many initial terms as possible are set to zero, the 
next term in the extremal weight enumerator is always positive (so the minimal 
distance of the hypothetical extremal code is precisely 4[n/24] + 4). 

Lengths that are multiples of 24 are especially interesting. Figure ^ shows 
the leading term in the extremal weight enumerator for a binary self-dual code 



of length n = 24m (this is sequence A34414). It is "obvious" that these numbers 



n 


Coefficient 





1 


24 


759 


48 


17296 


72 


249849 


96 


3217056 


120 


39703755 


144 


481008528 


168 


5776211364 


192 


69065734464 



Fig. 4. Number of codewords of minimal weight n/6 + 4 in extremal weight 
enumerator of length n. 



are growing rapidly, and in fact it is shown in J jj that if n = 24m then the 
leading coefficient is 

24m\ /5m - 2\ /4m + 4 
5 ) V m 5 

The next term in the extremal weight enumerator (the number of codewords of 
weight 4[ra/24] + 8) also grows rapidly — see Fig|| (sequence A34415). 



Again it is "obvious" that these numbers also grow rapidly, especially if 
one examines a more extensive table, where the number of digits continues to 
increase with each term. Yet — this shows the danger of drawing conclusions 
just from numerical data — it is proved in |6!| that this sequence goes negative 
at n = 3696 and stays negative forever. (Sec Theorem 29 of |82j for more precise 
information.) 
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n 


Coefficient 





1 


24 


2576 


48 


535095 


72 


18106704 


96 


369844880 


120 


6101289120 


144 


90184804281 


168 


1251098739072 


192 


16681003659936 



Fig. 5. Number of codewords of next-to-minimal weight n/6 + 8 in extremal 
weight enumerator of length n. 



So codes corresponding to these extremal weight enumerators certainly do 
not exist for n > 3696 (since the weight enumerator of a genuine code must 
have nonnegative coefficients). They are known to exist for n = 24 (the Golay 
code) and 48 (a quadratic residue code), but every case from 72 to 3672 is open. 
For further details, and in particular for a description of the analogous situation 
when n is a multiple of 8 but not of 24, see 

The situation for lattices is similar. Consider an even unimodular lattice A 
of dimension n. with thcta series 



Examples are the eight-dimensional root lattice Eg, with theta series 
/ = 1 + 240g 2 + 2160g 4 + 6720g 6 + • • • 

OO 

= 1 + 240 ^ct 3 (m)q 2m , 

m— 1 

where 173(771) is the sum of the cubes of the divisors of m, and the 24-dimensional 
Leech lattice, with theta series 

.9 = 1 + 196560g 4 + 16773120^ + 398034000<7 8 + • • • 

OO 

= f - 720<z 2 Y[ (1 q 2m ) 24 . 

m—l 

The Taylor series expansion of the last product defines the famous Ramanujan 



numbers. (The coefficients of / and g give seq uences [A4009| and [A8408; , respec- 
tively; the Ramanujan numbers form sequence A594j .) 



A theorem due essentially to Hecke says that the thcta scries of any such 
lattice A is a polynomial in / and g. We may then define extremal thcta se- 
ries just as we defined extremal weight enumerators. It is known [^39| that the 
leading term in the extremal theta series is positive, but that again the next 
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coefficient eventually becomes negative (even though the sequence begins with 
about 1700 exponentially growing terms). The dimension where the first nega- 
tive coefficient appears is around 41000. The leading terms are shown in Fig. ^| 
(sequence A34597). 



n Coefficient 

~~ r 

24 196560 

48 52416000 

72 6218175600 

96 565866362880 

120 45792819072000 



Fig. 6. Number of shortest vectors u (with u ■ u = n/12 + 2) in extremal thcta 
scries at dimension n. 



Lattices corresponding to the extremal theta series exist for n = 24 (the 
Leech lattice) and 48 (at least three inequivalent lattices exist), but every case 
from 72 to about 41000 is open. The 24- and 48-dimensional extremal lattices 
provide record kissing numbers in those dimensions; it would be nice to have a 
72-dimcnsional example. 

8 Levine's sequence 

In the summer of 1997 Lionel Levine [[32| submitted a new sequence to the table, 
a sequence of such beauty that it immediately caught the attention of several 
of my colleagues. It is constructed via the array in Fig. If a row of the array 



1 1 


1 2 


1 1 2 




1 


12 3 


111 


2 


2 3 4 


111122233 4 


4 


5 6 7 


111111122222233333 4 


1 




4455556667778899 10 11 12 13 14 



Fig. 7. The array that produces Levine's sequence 1, 2, 2, 3, 4, 7, 14, . . .. 
contains the numbers 

ai d2 0,3 . . . a,k-i a,k 
15 



then the next row contains 



a k l's, a fc _i 2's, a fe _ 2 3's, . . . 
Levine's sequence ( |A11784 ) is obtained by taking the last term in each row: 



1, 2, 2, 3, 4, 7, 14, 42, 213, 2837, 175450, 139759600, 6837625106787, , . 
266437144916648607844, 508009471379488821444261986503540, ... ' 

The terms grow unexpectedly rapidly! The nth term L n is 

(i) the sum of the elements in row n — 2 

(ii) the number of elements in row n — 1 

(iii) the last element in row n 

(iv) the number of l's in row n + 1 

Furthermore, if s(n, i) denotes the sum of the first i elements in row n, then we 
have 

(v) L n+2 = s(n,L n+1 ) 

L n+1 

(vi) L n+3 = ^2 s ( n '«) 

i=l 

(vu) L n+4 = 2^ 

i=l 



The latter identity was found by Allan Wilks 1 101 ] , who also found a more 
complicated formula for L n +5, and used it to compute the last two terms shown 
in (g). No other ter ms are known! 

As to the rate of growth, we have only a crude estimate. Bjorn Poonen and 
Eric Rains [MJ showed that 

logL n ~cr n , (3) 

where r = (1 + \/5)/2. Sketch of Proof, (a) L„+2 < L n +iL n , and so logL n is 
bounded above by a Fibonacci-like sequence, (b) The sum of the (n + l)st row 
is at most 



in 

2 



which implies 



2L n+ 2 2L„ + i 2L n 

and so log(L„+i/2L„) is bounded below by a Fibonacci-like sequence. 

Colin Mallows |i8| has determined numerically that a reasonably good ap- 
proximation to L n is given by 

— e C2T " 

Cl 

where c\ f» 0.277, c 2 ~ 0.05427. It would be nice to have better estimates for 
L n , and one or more additional terms. 
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9 Arrays of Numbers 



Besides number sequences, the database also contains arrays of numbers that 
have been converted to sequences. Triangular arrays are read by rows, in the 
obvious way. E.g. Pascal's triangle of binomial coefficients 

1 

1 1 
1 2 1 
13 3 1 
1 4 6 4 1 

becomes sequence A7318| : 

1,1,1,1,2,1,1,3,3,1,1,4,6,4,1,... 

Square arrays are read by antidiagonals, usually in this order: 

a a 2 a 5 a 9 . . . 
a\ a& a% . . . 
a 3 a 7 . . . 
ae . . . 



E.g. the Nim-addition table |TJ 



1 2 3 4.. 

1 3 2 5.. 

2 3 1 6.. 

3 2 1 7.. 



becomes sequence |A3987 : 

0,1,1,2,0,2,3,3,3,3,4,2,0,2,4,... 

Other classical arrays are the Stirling numbers of both kinds, Eulerian num- 
bers, etc. 

A less well-known array arises from Gilbreath's conjecture. This conjec- 
ture states that if one writes down the primes in a row, and underneath the 
absolute values of the differences, as in Fig. ||, then the leading terms (shown 
underlined ) of all rows except the first are equal to 1 (j33j, §A10). The corre- 
sponding sequence (A36262) is 



2,1,3,1,2,5,1,0,2,7,1,2,2,4,11,1,.. 
Odlyzko p6| has verified the conjecture out to 3 x 10 11 . 
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Fig. 8. Gilbraith's conjecture is that the leading terms of all rows in this array 
except the first are always 1 (the top row contains the primes, subsequent rows 
are the absolute values of the differences of the previous row). 



10 The Wythoff array 

This array shown in Fig. ^| has many wonderful properties, some of which arc 
mentioned here. I learned about most of these properties from John Conway 
|jl3| , but this array has a long history — see Fraenkcl and Kimberling [28], 
Kimbcrling ||,||,|g|,g|,|3,g8),||, Morrison @ and Stolarsky @, 
|)7| . It is r elated to a large number of sequences in the database (the main entry 
is [A35513D . 




1 
2 
3 

4 8 

5 9 

6 11 

7 12 



1 2 3 5 8 13 21 34 55 

4 7 11 18 29 47 76 .. 

6 10 16 26 42 68 .. 

9 15 24 39 63 .. 
12 20 32 52 
14 23 37 60 
17 28 45 73 
19 31 50 81 



Fig. 9. The Wythoff array. 



Construction (1). The two columns to the left of the vertical line consist respec- 
tively of the nonnegative integers n, and the lower Wythoff sequence ( |A20l| ), 
whose nth term is [(n + 1)t]. The rows arc then filled in by the Fibonacci rule 
that each term is the sum of the two previous terms. The entry n in the first 
column is the index of that row. 

Definition. The Zeckendorf expansion of a number n is obtained by repeatedly 
subtracting the largest possible Fibonacci number until nothing remains. E.g. 
100 = 89 + 8 + 3 = Fn + Fq + F4. The Fibonacci successor to n, Sn, say, is found 
by replacing each Fi in the Zeckendorf expansion by i^+i- E.g. the Fibonacci 
successor to 100 is S100 = F u + F 7 + F 5 = 144 + 13 + 5 = 162. 
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Construction (2). The two columns to the left of the vertical line in Fig. ||read 
n, 1 + Sn; then after the vertical line the row continues 



to Sm SSm SSSm SSSSm 



where m — n + 1 + Sn. 



Construction (3). Let { 51, 52, S 3, . . .} = {2, 3, 5, 7, 8, 10, 11, . . .} be the sequence 
of Fibonacci successors (A22342). The first column to the right of the line consists 
of the numbers not in that sequence: 1, 4, 6, 9, 12, . . . ( A7067). The rest of each 
row is filled in by repeatedly applying S. 



Construction (4). The entry in row n and column k is 

[(n + l)rF fe+2 ] + F k+1 n 

(where k = indicates the first column to the right of the vertical line) . 
Some properties of the array to the right of the line are the following: 

(i) The first row consists of the Fibonacci sequence 1, 2, 3,5,8,... 

(ii) Every row satisfies the Fibonacci recurrence. 

(iii) The leading term in each row is the smallest number not found in any earlier 
row. 

(iv) Every positive integer appears exactly once. 

(v) The terms in any row or column are monotonically increasing. 

(vi) Every positive Fibonacci-type sequence (i.e. satisfying a(n) = a(n — 1) + 
a(n — 2) and eventually positive) appears as some row of the array. 

(vii) The terms in any two adjacent rows alternate. 

There are infinitely many arrays with Properties 1-7, see Ej . 
The nth term of the vertical para- Fibonacci sequence 

0, 0, 0, 1, 0, 2, 1, 0, 3, 2, 1, 4, 0, 5, 3, 2, 6, 1, 7, 4, 0, 8, 5, . . . 

( IA19586D gives the index (or parameter) of the row of the Wythoff array that 
contains n. This sequence also has some nice fractal-like properties: 

(a) If you delete the first occurrence of each number, the sequence is un- 
changed. Thus if we delete the underlined numbers from 

0,0,0,1,0,2,1,0,3,2,1,4,0,5,3,2,6,1,7,4,0,8,5,... 

we get 

0, 0, 0, 1, 0, 2, 1, 0, 3, 2, 1, 4, 0, 5, 3, 2, 6, 1, 7, 4, 0, 8, 5, . . . 

again! 

(b) Between any two consecutive 0's we see a permutation of the first few 
positive integers, and these nest, so the sequence can be rewritten as (read across 
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the rows): 




1 

2 1 

3 2 1 4 

5 3 2 6 1 74 

085392 10 61 11 74 12 

The nth term of the horizontal para- Fibonacci sequence 

1,2,3,1,4,1,2,5,1,2,3,1,6,1,2,3,1,4,1,2,7,1,2,. 



( IA35612D gives the index (or parameter) of the column of the Wythoff array that 



contains n. This sequence also has some nice properties. 

I hope I have said enough to convince you that the Wythoff array is well 
worth studying and full of surprises. 



11 The Boustrophedon transform 

The Taylor series for sin a; and cos a; are easily remembered, but most people 
have trouble with 



tana; = Ix + 2— + 16— + 272— + • • • , 

2 4 6 

x x x 
secx = 1 + l_ + 5 ¥ +61- + --- . 

However, their coefficients can be calculated from the array in Fig. nu. The 
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Fig. 10. The secant-tangent triangle. 



nonzero entries on the left arc the secant numbers (|A364| ) 



1, 1, 5, 61, 1385, 50521, 2702765, 
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and those on the right are the tangent numbers (A182) 



1, 2, 16, 272, 7936, 353792, 22368256, . 



while the combination of the two sequences (Alll): 



1, 1, 1, 2, 5, 16, 61, 272, 1385, 7936, ... (4) 

are usually called the Entringer numbers. The latter count permutations of 
{1,2, ... ,n} that alternately fall and rise. 

This array is filled by a rule somewhat similar to that for Pascal's triangle: 
the rows are scanned alternately from right to left and left to right, the leading 
entry in each row is 0, and every subsequent entry is the sum of the previous 
entry in the same row and the entry above it in the previous row. (This is the 
boustrophedon or "ox-plowing" rule.) The earliest reference I have seen to this 
triangle is Arnold who calls it the Euler-Bernoulli triangle. However, it may 
well be much older origin. |f72| gives many other references. 

Richard Guy |m| observed that if the entries at the beginnings of the rows 
are changed from 1, 0, 0, . . . to say 1, 1, 1,1,1,... or to 1, 2, 4, 8, 16, . . . then the 
numbers that appear at the ends of the rows form interesting-looking sequences 
that were not to be found in |k|, and asked if they had a combinatorial inter- 
pretation. Using 1,1,1,... for example the triangle becomes 

1 

1 2 
4 3 1 
15 8 9 
24 23 18 10 1 
1 25 48 66 76 77 



yielding the sequence (A667) 



1,2,4,9,24,77,294,1309,... (5) 

We may regard this process as carrying out a transformation (the boustrophe- 
don transform) of sequences: if the numbers at the beginnings of the rows are 
ao, <xi, ci2, . . • (the input sequence) then the numbers at the ends of the rows, 
bo, b\, &2, ■ ■ • (say) arc the output sequence. In j] | we showed that there is a 
simple relationship between the input and output sequences: their exponential 
generating functions 

00 x 11 °° x n 

•A{x) = Va„- r, 8(i) = Vi)„- 
z — < n \ £ — t n \ 

n=0 n=0 

are related by 

B{x) = (sec x + ta,nx)A(x) . 
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We also give a combinatorial interpretation of the {&„}. E.g. in (||), b n is the 
number of up-down subsequences of {1, ... ,n}, so that 63 = 9 corresponds to 
0,1,2,3,12,13,23,132,231. 

The Entringer sequence ([lj) then has the property that it shifts one place 
left under the boustrophedon transform. The le xicogra phically earliest sequence 
that shifts two places left under this transform ( A661 ) is 

1, 0, 1, 1, 2, 6, 17, 62, 259, 1230, 6592, . . . 

We do not know what this enumerates! 

Many examples of similar "eigen-sequences" for other transformations of se- 
quences can be found in Donaghey |23|, Cameron ||, and espe cially B . E.g. the 
sequence giving the number of planted achiral trees ||(J , | | ( A323§| ): 

1,1,2,3,5,6,10,11,16,19,26,... 

has the property that it shifts left one place under the "inverse Mobius trans- 
formation" given by 

b„ = y^«d . 

iiin 



12 Tchoukaillon solitaire (or Mancala, or Kalahari) 

These are ancient board games, with hundreds of variants and many different 
names. The version to be described here is called Tchoukaillon solitaire. It has 
been studied by several authors (see for example Betten Q and Brolinc and 
Loc b pi ). It is played on a board with a row of holes numbered 0, 1, 2, . . . (see 
Fig. In]) . 



5 4 3 2 1 

Fig. 11. A move in Tchoukaillon solitaire. 
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The game begins with n stones placed anywhere except in hole 0. A move 
consists in picking up the stones in some hole and placing one in each lower- 
numbered hole. If the last stone falls in hole then play continues, otherwise the 
game is lost. The objective is to get all the stones into hole 0. 

The game is interesting because there is a unique winning position for any 
number of stones. These winning positions are shown in Fig. [l^ (sequence A28932), 
and can be found by playing the game backwards. 



n 


Position 








1 


1 


2 


20 


3 


21 


4 


310 


5 


311 


6 


4200 


7 


4201 


8 


4220 


9 


4221 


10 


53110 


11 


53111 


12 


642000 


13 


642001 



Fig. 12. The unique winning position for n stones in Tchoukaillon solitaire. 



The array can be more explicitly constructed by the rule that if the first in 
a row (counting from the right) is in position i, then the next row is obtained by 
writing i in position i and subtracting 1 from all earlier positions. The sequence 
of successive values of i (A2892C) is 



1,2,1,3,1,4,1,2,1,5,1,6,1,2,.. 



Let t(k) denote the position where k occurs for the first time in this sequence. 
The values of i(2), i(3), ... are (sequence A2491): 

1,2,4,6,10,12,18,22,30,34,42,... 

This sequence has some very nice properties. It has been investigated by (in 
addition to the references mentioned above) David [jlTj, Erdos and Jabotinsky 
[I and Smarandache p4j . 
(i) t(n) can be obtained by starting with n and successively rounding up to 
the next multiple of n — 1, n — 2, . . . , 2, 1. E.g. if n = 10, we obtain 

10 -> 18 -> 24 -> 28 -> 30 -> 30 -> 32 -> 33 -> 34 -> 34 , 
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so rj(10) = 34. 

(ii) The sequence can be obtained by a sieving process: write 1,2, .. . in a 
column. To get the second column, cross off 1 and every second number. To 
get the third column, cross off the first and every third number. Then cross off 
the first and every fourth number, and so on (see Fig. [l3|) . The top number in 
column n is t(n). Comparison of Figures |l2| and [l3] shows that connection with 
the solitaire game. 
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Fig. 13. A sieve to generate the sequence i(l), rj(2), . . . = 1, 2, 4, 6, 10, 12, 18, . . .. 
At stage n, the first number and every nth are crossed off. 



(iii) Finally, Broline and Loeb |6| (extending the work of the other authors 
mentioned) show that, for large n, 

o 

tin) = h 0(n) . 

It is a pleasant surprise to see 7r emerge from such a simple game. 
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